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An often ignored aspect of electromagnetic radiation from 
antennas is the characterisation of their near-fields. A 
computer program, Numerical Electromagnetics Code (NEC), is 
validated for accurate near-field computations and applied to 
a model of a broadcast monopole. E- and H-fields are plotxed 
as a function of position along the antenna for various 
distances from the surface. The fields are also plotted as 
a function of radial distance outward for various heights. 
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INTRODUCTION 



A. PROBLEM ENVIRONMENT 

Maxwell's equations describe, for all space, the Electric 
(E) and Magnetic (H) Fields that are the riedium for trans- 
mitting information (vix'tually) instantaneously/ over 
extended distances. This practical application of electro- 
magnetic (EM) radiation has generated considerable interest 
in its far-field* characteristics. The E and H vectors can 
be analytically determined when the current (I) on the 
radiator (from which the field vectors are generated) has 
been described. For the far-field region of monopoles and 
dipoles, acceptable results are obtained by assuming a 
sinusoidally varying current. However, when the field 
vectors are calculated for field points at distances less 
than \/2tt (the near-field region), significant errors occur 
between calculated values and empirically measured values . 
This discretancy is not she result of inaccuracies in 
Maxwell’s Equations, but can be wholly attributed to an 
inaccurate description of the current. 

*The far-field, for a 7r*oncpole or dipole, is defined as dis- 
tances from the radiation greater than a/2t, where \ is th 
wavelength. For free space, A=c/f where f=the frequency a 
which the radiator is excited, and c=the speed of lighr in 
free space. 
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In the past, the inaccuracies of the near field calcula- 
tions were not a problem — there was no practical application 
to which the energy could be applied. However, in recent 
years there has been an escalating interest in the close-in 
problem. This interest has been generated by an increasing 
proliferation of devices which use Radio Frequency (RF) energy 
in such a manner that the user is exposed to near-field 
radiation (e.g. microwave ovens, hand-held walky-talkies , 
wireless telephones, RF sealers, etc.). 

The U.S. Navy was one of the first organizations to 
express an interest in these fields. The interest was the 
result of a concern for the hazards these fields might pose 
to fuels, ordnance, other EH operational and test equipment, 
and to personnel. The hazards posed by fuels, ordnance, 
and other EM equipment result in overt effects; the hazards 
posed to personnel are not so obvious. 

With virtually unanimous consensus, the literature 
identifies the heating of body tissue as the primary hazard 
to personnel. However, some distinctly non-thermal effects 
have been observed. Some of these effects [Ref. 1] include: 

a) minor changes in blood properties 

b) a "buzz” heard by certain people when exposed to 
microwave radiation* 

*The "buzz" is assum^ed to be a function of the pulse repetition 
rate (?RR) rather tnan the carrier frequency. 
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c) abnormalities of the chromosome structure 

d) movement, orientation, and polarization of protein 
molecules 

e) changes in the transport rate at the blood-brain 
barrier 

f) comfort imbalances, such as epigastric distress, 
emotional upset, and nausea. 

The significance of these effects is not yet understood. 

The general public has not been exposed to near-field 
radiation hazards of sufficient intensity or duration to be 
cause for alarm; however, there are certain vocations in 
which the incumbents have an elevated probability of 
exposure to this type of field.* 

B. THESIS STATEMENT, SCOPE AND LIMITATIONS 

In order to provide more accurate data for people 
researching biological effects of near-field radiation, this 
study will validate a computer program which computes the 
near electric, miagnetic, and composite (peak) fields of an 
arbitrary radiator. The program v;ill also be used to 
compute these near fields for a model to be described below. 
The progr'am. Numerical Electromagnetics Code (NEC) , was 
developed at Lawrence Livermore Laboratory under the 



'•Sailors a~ sea or steeplejacks who effect antenna repairs , 
for instance. 
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sponsorship of the Naval Ocean Systems , Center (NCSC) and the 
Air Force Weapons Laboratory. It can be easily learned and 
applied directly by a researcher who has some knowledge of 
electromagnetic theory. 

The model to be used for xhe investigation will be a 75 
meter high broadcast monopole radiating above a perfectly 
conducting ground plane. The ground plane is located in the 
X-Y plane; the monopole is co-axial with the Z-axis . The 
antenna will be excited at its base with a volxage sufficient 
to produce a radiated power of 1000 watts. The excitation is 
at a frequency of 1 Megahertz, which has a wavelength of 300 
meters; hence, the monopole is a quarter-v/avelength radiator 
operating ar resonance. A diagram of the model is shown in 
Figure 1. 

In simulating this model, NEC requires that the antenna 
be broken into virtual short straight segmenrs. In 
consonance with this requirement, three different configura- 
tions will be used for computations; a model composed of 5, 
15, and 25 segmenrs of 15, 5, and 3 meters, respectively. 

Using cylindrical coordinates, the radial and parallel 
electric fields, and the rangential magnetic field will be 
computed. Field values will be computed in .1 merer 
increments for points parallel to the z-axis at distances of 
1", 2, 5, 10, 100, and 1000 radii (1 radius = .3 meters). 



*1 radius = surface of the moncpole. 
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Field values will also be calculated along the radial axis 
(Y-axis) in increments of .5 meters for heights (h) of 2h, 
Ih, l/2h, l/4h, and 2/lOOths h.-- 

The near-field examined in this study is the cylindrical 
volume bounded by a cylinder of one wavelength radius 
centered on the Z-axis , rising from the ground plane to a 
height twice that of the monopole. 



This is the approximaxe height of an average mans upper 
torso . 
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II. DESCRIPTION OF SELECTED ANALYSES 



A. CLASSICAL ANALYSIS 

The antenna utilized in this study is a member of the 
class identified as *thin-wire’ antennas. The criteria 
which adumbrate the thin-wire approximation are [Ref. 2]: 

a) Transverse currents are negligible relative to axial 
currents 

b) Circumferential variations of the axial current are 
negligible 

c) The current can be represented by a filament current 
on the wire axis 

d) Boundary conditions on the electric field need be 
enforced only in the axial direction. 

By far, the largest subset of thin-wire antennas is a 
class known as short antennas. Short antennas are defined 
as having one dimension (h) much greater than the other 
dimensions, but h is also much less than a wavelength, 
typically h < \/8. Because of the predominance of this 
sub-set, it is reasonable that it should be used as a model 
for classical analysis. 

Based on Maxwells Equations and geometric considerations 
as shown in Figure 2, the following non-zero complex scalar 
field components can be derived: 
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( 1 ) 
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^ — cos 0(0 + -^- 

r, 




( 2 ) 




( 3 ) 



where; n = rhe intrinsic impedance of the medium, I = antenna 
currenx , \ - wavelength, S = 2ir/X = phase constant, h = 
monopole height, r = distance of field observation point 
from monopole base. Since there is no variation in $: 



The interpretation of these equations under classical 
analysis require two additional assumptions: 

a) The distance to the field observation point is large 
relative to the moncpole height, and 

b) An harmonically time varying current is uniform over 
the height of the monopole. 

The first of these assumptions is made to focus attention to 
that volume of space where the lines of force have detached 
themselves from rhe anrenna and are approaching the form of 
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a plane wave front. The second assumption is to form a 
first approximation of the current distribution on the wire 
which is solvable in closed form. Extensive analysis over 
many years has shown this to be a valid assumption when the 
antenna is excited at a single point. 

In this study, we will be interested in fields where the 
observation points are not at distances considered large 
relative to the monopole height. When This is the case, the 
diagram of Figure 2 must be modified as shown in Figure 3*. 
From the geometry of this figure, the following non-zero 
complex scalar field components can be derived: 
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^Equations (1) through (3) are expressed in spherical coor- 
dinates while Equations (5) through (7) have been Transformed 
TO cylindrical coordinates. 
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A current of the form: 



I(z) = I^sinSh for 0 <_ z ^ h (8) 

has been inserted into Equations (5) through (7). 

Equations (5) through (7) indicare that each field 
component is the phasor sum of contributions from three point 
sources. The first source is located at the top of the mono- 
pole, a second image source is located in a symmetric position 
when the first source is reflected through the ground plane. 
The third source is located at the base of the monopole; its 
relative amplitude is modified by the (-2cos6h) factor. 

Since we have a solution for computing the field points, 
it is now reasonable to determine the range over which it is 
valid. By inspection we note that for r = \/2 it in Equations 
(1) through (3), the magnitude of each non-zero rerm in the 
parenthesis is equal. This implies that r = X/2-:i may be a 

transition point. Further investigation indicaxes that for: 

3 . 2 

a) r>,\/2T7: the 1/r terms in E„ and H. , and the 1/r 

rj (L 

term in E are dominant 
r 

3 2 

b) r<X/27r: the 1/r terms in E. and E , and the 1/r 

9 r 

term in H. are dominant. 

Since we know that power drops off as the inverse square of 
distance, and that the power flow at a field point can be 
defined in terms of its E-field as: 



(9) 



P = E^/n 

2 2 

we can extrapolate from a) above (E .8/n« k/r . where k = a 
constant) that the far-field begins at r = a/2tt and continues 
out to infinity, while the near-field begins at the monopole 
and continues out to r = X/2 tt.* 

It is interesting to note the complementary relationship 
that is maintained by the power in this transition region. 
Instantaneous power flow at a point in space is given by the 
complex Poynting vector: 

S* = 1/2(E X H) (10) 



For the short monopole, the complex Poynting vector is 
oriented in the radial direction and given by: 




n (Ih 
2 4tt 



)(sin9)e”^^^)^(^ 




( 11 ) 



The real part of this equation gives the average rate of flow 

2 

of Real power («l/r ) while the Imaginary part is proportional 
to the difference between average magnetic energy density 



"This argument is offered not as a proof, but as a reasonable 
rationale which can be empirically/’ verified. 
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and average electric energy density («l/r ). Interpretation 
of this equation indicates that for: 

a) r>\/2ir: Reactive power flow decreases rapidly 

1 ^ 

Real power flow decreases not so rapidly 

r 

b) r = X/2it: Average Real power flow and average Reactive 
power flow are equal 

c) r<X/2-iT: Reactive power flow becomes much greater 

5 2 

(*r ) than Real power flow (®r ). 

The transverse component of the complex Poynting vector: 



S- 



rlh ^-3Sr . 2 . . . ,S 

= nCri — e ) cos9s me ( -i (-^ 
4 77 O 



■) ) 9 



( 12 ) 
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is completely Imaginary, indicating no Real power flow, and 

3 

it falls off rapidly with increasing distance due to 1/r and 
1/r^ contributions.* 

The deviant behavior of the fields encountered at 
distances r<X/27r is the result of moving from a region vjhere 
the lines of force are closed on themselves and approaching a 
plane wave front, to a region where the lines of force still 
maintain a strongly splaerical wave front and some of the lines 
are still attached to the antenna. Even in the simplest of 

*Also noted is the fact that the Electric field changes from 
Elliptical polarization to Linear polarization in passing 
from r>X/ 277 to value of r<X/27r when 9 0° , 90°, or 18 0°. 

For these angles the polarization of the electric field is 
linear on both sides of r = X/2 t7. 



cases, the current distribution in conductors is complex. 

This complexity is the result of non-uniform currents and 
charge densities, capacitance and inductance generated by 
these non-uniformities , and discontinuities from impressed 
forcing functions. 

It is well established that as the fi.eld observation point 
moves within the x/2tt distance from the antenna, that the 
measured values deviate from those predicted by Equations (5) 
through (7). It is also well established that Maxwells 
equations are valid for all space. Consequently, the con- 
clusion is That the sinusoidal current distribution is not 
sufficiently accurate for determining field poinrs within 
the near-field region. The problem then is to find a better 
description of the current distribution. This is the domain 
of contemporary analysis. 

Before proceeding, it is necessary to reflect thar the 
analysis performed thus far has been for a short antenna. 

The model for this study is a resonant monopole; when the 
antenna is lengthened to resonance some interesting 
characteristics appear. Primary among these is the fact 
thar the reactive ccm.pcnent of the impedance vanishes , which 
says the input reactance goes to zero. The current on the 
antenna and the voltage required to place it there are 
determined solely by the resistance, the real part of the 
impedance. When a sinusoidal current is assumed on the 
antenna, the tangential field component is given by: 
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(13) 



The (-2cosSh) term representing radiation from the base of 
the monopole has disappeared. The monopole appears as a 
point source radiating from its height, and an image source 
found at the reflection through the ground plane. 

Also, the resonant condition occurs at a length less 
than that expected from the X/4 relationship. This is the 
result of an "effective" lengthening of the anxenna by 
capacitance and "fringing" at the end-cap of the wire. 
Capacitance is caused by a build-up of charge at the end of 
a wire of finite width. The change in direction of the 
geometry at the end of the wire also changes the direction 
of the lines of force emanating perpendicular from it; this 
is the fringing effect. These factors cause resonance to 
occur for an antenna length ofh=X/4-o where 6 is a 
function of the capacitance caused by the two effects 
discussed above. 

The short monopole development discussed above is very 
restrictive and cannot be applied directly xo resonant 
monopoles. However, when the resonant peculiar characteris- 
tics mentioned are taken into consideration, the resulting 
development will follow the same orocedure as used for the 
short monopole. For the resonant moncpole development and 
results, see Ref. 4. 
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CONTEMPORARY ANALYSIS 



Schelkunoff [Ref. 3: pp. 370-374] identifies three methods 
for obtaining a closer approximation to the actual current 
distribution on an antenna. They are: 

a) The Integral Equation method 

b) The Sweep-Off method 

c) The Mode Theory of Antennas method. 

Each are discussed below, in succession. 

1 . The Integral Equation Method 



The assumption is made that the currenx is distributed 



on the surface of a hollow cylinder which is divided into 
filaments of angular density (I(z*)/2t:) and angular width 
(d<5'). The Kernel Function is modified to allow for 
variations in and is expressed as an integral around the 
circumference of the wire as : 




(14) 



where : ^ 

4irr 



-i 3r 

0 



wave runcTion. 



The field intensity for a straight current filament: 




( 15 ) 



is modified by inclusion of Equation (14), with the boundary 

condition: (E + E^ = 0) enforced to give: 

z z 



-E^(z) 

2 





d Z 



+S^^) I(z') dz' 



( 16 ) 



This Integral Equation (from which the name of xhe method 
derives) is the limit (as ) of a series of equations of 
the form: 

+ 6^ (zz)] Az I(z ) 
n ’ n n 

(17) 



-E^(z) = 
z 



n = l - 



( z 



n 



3z 



which is also of the form: 

-E^(z) = Z I(z) 
z mn 



(13) 



Equation (16) is a ’circuit equation' for an antenna and can 
be solved, in a manner analogous to the use of Kirchhoffs 
Equations for lumped networks, to an arbitrary degree of 
accuracy through the use of Equation (17). 
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It is instructive at this point to pause and 
consider how accurately our initial approximation of the current 
approached reality. Integrating Equation (IS) by parts and 
setting -E^(z)=0, as it is on the surface of the antenna, we 
obtain: 



(I(z')^(z,z’ ) 



+ I(z' ) 




z 2 

- ; ^ (^-^ + 3^I)^(z,z’) dz’=0 

z^ dz ' 



( 19) 



This equation shows that as the antenna radius approaches 

zero, the first term of Equation (19) approaches a constant 

limit for all z, except z = or z = Z2* Sufficiently far 

from the ends this constant is sm.all. The ^(z,z’) factor in 

the second term is infinite at z = z ’ and large in its 

2 

cl ^ 2 

vicinity; hence, the ( — ^ + S I) factor must be small. 

dz ' " 

Therefore, as the radius of the antenna approaches zero, 
d^I 2 

( y+S I) approaches zero, which implies I(z) approaches 

dz ’ 

a sinusoidal form. 

2 . The Sweep-Off Method 

This method starts with Equation (15), the tangential 
field equation for a filament current, being integrated twice 
by parts. Then an approximation to the current distribution 
on a transmitting antenna of finite radius, assum.ed to be 
given by: 
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I(z ’ ) 



A cosBz ' + B sinSz ' 



( 20 ) 



is substiruted into the integrated form of Equation (15) to 
yield: 



-jBr, 



E, = [I’(z.)^ 

z 4iT(DE 2 Vr 



-jSr^ -jSr^ 



) 



-j8r^ 



( 21 ) 



This equation is an exact expression for the electric field 
intensity parallel to a sinusoidal current filament. The 
antenna is then subjected to a compensating impressed field, 
-E^ , and the correction terms from the first aDoroximation 
are computed. The name for this method results from the 
concept of ’sweeping-off’ the surface of the antenna the 
residual tangential electric intensity by applying an equal 
and opposite intensity. 

It should be noted here that alrhough Equation (21) is an 
exact expression for E , that E does not vanish along the 
filament. This apparent dichotomy results from the erroneous 
assumption of uniform convergence for E^. In actuality, the 
convergence of E is non-uniform at p = 0 , where p is the 
distance from the filament. For an antenna radius r = a. 
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there exists a remainder term 1 ^( 2 ) in addition to I(z); 

It(z) is sufficient to cause E (a) =0. As a^O , I,(z)->0: 

1 z 1 

consequently, its contribution to E^ for any fixed value of 
p, greater than a, will approach zero with a. In the limit, 
E^ is given by Equation (21) for p>0, but at p = 0 Equation 
(21) does not hold. 

3 . Mode Theory of Antennas 

This method is based on the solution of i!axwells 
Equations subject to boundary conditions at the antenna 
surface, and the surface of the source of power. The method 
consists of calculating modes of propagation consistent with 
the boundary conditions at the lateral surfaces of the 
antenna, then the modes consistent with free-space propaga- 
tion. These modes are combined to satisfy all boundary 
conditions and the current is obtained as the sum of two 
components: the TEM or Principal wave, and the Complementar'y 

wave consisting of all higher order modes. The principal 
mode is governed by ordinary transmission line equations 
(in terms of distributed capacitance and inductance). 

4 . Discuss ion 

Of the methods discussed above, 'Sweeping-Off’ is the 
most elementary and consequently, the least precise. The 
Mode Theory of Antennas excites a current in the conductor 
for the TEM mode, and a current for the complementary wave, 
which is a sum of currents, one for each of the higher order 
modes excited. The excited modes are strongly dependent on 
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